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Abstract. We describe here some new results concerning the Lorentzian Barrett-Crane model, a well- 
known spin foam formulation of quantum gravity. Generalizing an existing finiteness result, we provide 
a concise proof of finiteness of the partition function associated to all non-degenerate triangulations of 4- 
manifolds and for a class of degenerate triangulations not previously shown. This is accomplished by a 
suitable re-factoring and re-ordering of integration, through which a large set of variables can be eliminated. 
■ The resulting formulation can be interpreted as a "dual variables" model that uses hyperboloid variables 

associated to spin foam edges in place of representation variables associated to faces. We outline how this 
method may also be useful for numerical computations, which have so far proven to be very challenging for 
£SJ . Lorentzian spin foam models. 

m ■ 

1. Introduction 

Spin foam models offer a promising new quantum mechanical picture of space-time geometry. The spin 
foam program has been developed as a path integral formulation of loop quantum gravity 1 . Of the spin foam 
models proposed in the literature, perhaps the most studied have been the Barrett-Crane models. Originally 
constructed using representations of the group Spin(A) associated to manifolds of Riemannian signature 0, 
a subsequent model given by Barrett and Crane in [§] based upon representations of the Lorentz group 
SO(3, 1) is considered to be more physically realistic. 

While important strides have been made in numerical computations within the Riemannian framework 
0111 E]> numerical computations with the Lorentzian model have proven much more difficult. An important 
difference between the models lies in the space of representations, which in the Lorentzian case is continuous 
rather than discrete. In addition, the homogeneous space used to define vertex amplitudes is non-compact 
q-( in the Lorentzian case. 

The main result of the present work is that the Lorentzian partition function can, by a process we shall 
(3JT), refer to as "face factoring", be rearranged so that integration over the representation variables is performed 

exactly. After applying this transformation, computational or analytic effort can focus entirely on the 
homogeneous space integrals and the one-dimensional integrals associated to each edge of the spin foam. 

We note that the model that results from face factoring is strongly in the spirit of the dual variables 
picture proposed by Pfeiffer in {2Qj. While (20] dealt specifically with a Riemannian version of the Barrett- 
Crane model, our face-factored model can be understood as the realization of a dual variables picture in the 
Lorentzian case. 

The organization of the paper is as follows. In the next section, we recall the form of the Lorentzian 
Barrett-Crane amplitude and describe the face factoring method. Included in this section is an illustration 
of the method applied to the partition function of a simple 2-complex. In Section 3, we prove that the 
partition function associated to triangulations of a certain type is absolutely convergent in the face factored 
formulation, which justifies the interchange of integration required by the method. It follows that the 
original partition function is finite for triangulations of this type, which includes both the non-degenerate 
triangulations proven finite in Jlj as well as certain degenerate triangulations that have not previously been 
proven finite. It should be noted here that in the present work and in ^3], finiteness is shown for a specific 
choice of edge and face amplitudes; throughout this work our main interest is in the choice 2 due to Perez and 
Rovelli QHl- We also show how the method can be applied to an explicitly causal version of the Barrett-Crane 
model due to Livine and Oriti for certain types of triangulations. In Section 4, we briefly discuss the 

1 The reader is referred to [2] for a recent review of loop quantum gravity, and 1171 for a review of the spin foam program 
specifically. 

2 The physical content of these amplitude choices remains controversial; for example, see [§] for a critical discussion of this 
issue. In Section 3, we examine how certain changes in the face amplitude affect the type of degenerate triangulations for which 
we can show finiteness. 
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potential for numerical applications of the method and end with some concluding remarks. A bound required 
in Section 3 is proven in Appendix A, and in Appendix B we show that the finiteness results of Section 3 (in 
which a closed 4-manifold is assumed for clarity) also hold in the case of a 4-manifold with boundary. 

2. Description of the Face Factoring Method 

Review of the Lorentzian Barrett-Crane Model. To begin our discussion, we recall the form of a 
Lorentzian spin foam model. Spin foam models are defined on simplicial 2-complexes constructed from 
a set of vertices V, a set of edges E, and a set of polygonal faces F. Because our present work is with 
the Barrett-Crane model, we restrict ourselves to those 2-complexes which are the dual 2-skeleton of a 4- 
manifold triangulation 3 . The dual 2-skeleton is formed as follows: for every 4-simplex there is an associated 
dual vertex, and any two dual vertices whose associated 4-simplices meet at a tetrahedron are connected by 
a dual edge. For every face of the triangulation, the closed loop of 4-simplices that share that face gives rise 
to a two-dimensional polygon in the dual 2-skeleton, the vertices of which are dual to the 4-simplices and the 
edges of which are dual to the tetrahedra that are shared between neighboring 4-simplices. An important 
property of this construction is that the dual polygonal faces are in one-to-one correspondence with the faces 
of the original 4-manifold triangulation. Further background, including illustrations of dual skeleta, can be 
found in 0. 

We now introduce the Lorentzian Barrett-Crane model, which assigns to each dual 2-complex A a partition 
function Za as follows: 

(2.i) z A = r ■ ■ ■ r I n A f ) ( n a* ) ( n <m n $ 

> < Jo , \feF J Vee-B / Wv I feF 
feF 

where factors of the form pj arise from the measure on the principal series of SO(3, 1) representations; 
the amplitudes Af, A e , and A v will be defined shortly. For a given 2-complex, the state space associated 
with this partition function is the product space of representation variables pf (one for each face) and the 
multiple integration over this space is thought of as a path integration over all spin configurations. To 
find the transition amplitude between two spin networks (the basis states of 3-geometry in loop quantum 
gravity), a sum of Za over all 2-complexes interpolating between the spin networks is required. In the 
sum over 2-complexes, the amplitude Za for each 2-complex will be weighted by an appropriate measure to 
provide a normalization and possibly to regulate which 2-complexes can contribute. The choice of measure 
and questions of convergence have yet to be resolved — progress is likely to depend on what results can be 
found for the partition function of individual 2-complexes, which is the focus of our present work. 

In the expressions to follow, we write J| to indicate a formal product of symbols such as integral signs or 
measures. When a relation such as / 3 v appears below a product symbol (multiplicative or formal), the 
product is taken over all of the objects on the left hand side of relation that satisfy the relation; for example, 
in the f 3 v case the product would be over all faces / such that v is a member of /. The symbols v, e, and 
/ will always denote members of the sets V, E, and F, respectively. 

We now turn to the definition of the amplitudes A v , A e , and Af that appear in (2.1). Up to a regularization 
that will be defined below, the Lorentzian Barrett-Crane vertex amplitude of [Hj is defined as 



\e3v jH + / fBv 



where the kernel function K PJ is 



(2.3) K Pf (x,y) = 

The kernels K Pf are functions of the variables x e _ (f, v ), ^ e+ (/,u) £ H+ associated to the two dual edges that 
both contain v and are contained in /; the + and — serve to distinguish the two tetrahedra within a 4-simplex 
that share a face, but otherwise have no significance — any fixed convention can be chosen. As there are 
ten faces in the 4-simplex to which a vertex is dual, the integrand is a product of ten such kernels. The 



3 In Section 3, we extend the strict definition of triangulation to include degenerate triangulations which can be associated 
to topological 4-manifolds. 
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kernel function K Pf depends on the hyperboloid variables through their hyperbolic distance <j> on H+ , which 
is defined as 

(2.4) <P(x,y) = cosh~ 1 (x-y) 

for x,y G H+, where = {x e M 4 | x ■ x — 1, xo > 0}; here x ■ y is the Minkowski inner product 
x-y = x y a - x\y\ - x 2 y2 - x 3 y 3 . 

Although the expression for the vertex amplitude given above is generally infinite, we shall regularize it 
by fixing the value of one variable and dropping the corresponding integral over — the answer is 
independent of the choices [6J. 

It should be observed that the future 3-hyperboloid H\ is a homogeneous space of SL(2, C) and, unlike 
the homogeneous space S 3 that arises in the Riemannian model, it is non-compact. The other basic difference 
with the Riemannian model is that the relevant representations 4 of SL(2, C) are indexed by a continuous 
parameter pf, in place of a summation for each face the Lorentzian model requires an integral. 

While the original model of jg] specifies the vertex amplitude A v as given in (2.2), it leaves unspecified the 
edge and face amplitudes A e and Af. A proposal introduced by Perez and Rovelli in [IBJ and considered in 
further work chooses 5 Af — 1 and A e — 64 (pi, P2 1 Pz 1 Pi) for edges in the interior of the 2-complex. 

The function 64(pi,P2,P3,P4) is known as the eye diagram and can be defined as follows: 

,„ ^ , v 2 f°° sm(p 1 r e )sm(p 2 r e )sm(p 3 r e )sin(p 4 r e ) 

(2-5) e 4 (pi,P2,P3,P4j = / . ; 2 . , dr e , 

KP1P2P3P4 Jo smh (r e ) 

where Pi, • • • ,P4 are the spin variables labelling the four faces containing e. 

The results we have obtained with the face factoring technique have primarily been for the amplitudes 
Af — 1 and A e = 8 4 (pi)P2)P3)P4); we shall henceforth refer to this choice as the Perez-Rovelli Model. 

The Face Factoring Method. The major difficulty in computing Za as presented in (2.1) relates to the 
computation of the vertex amplitude (2.2). For each choice of the tenp/ variables, calculating a single vertex 
amplitude requires extensive computational effort in the form of Monte Carlo or quasi Monte Carlo integra- 
tion. In the absence of an exact expression or some other significantly more efficient means of calculation, it is 
prohibitively expensive computationally to evaluate the partition function for even the simplest 2-complexes. 
This situation has provided the main motivation for finding an alternative approach to calculating Za, in 
which no vertex amplitudes are explicitly computed. 

The idea of the present work is to perform the integration over the spin variables pf first, leaving a 
formulation in which the states to be integrated over are the configurations of H+ variables originating with 
the vertex amplitude integrand and the r e variables originating with the edge amplitude integrand. This 
reverses the conventional integration order implicit in (2.1) that has been adopted in most computational 
approaches to date. For example, in computations for the Riemannian case [H] the vertex amplitudes (equal 
to integrals over the homogeneous space) are computed for each spin configuration and the results summed 
to give the partition function. 

Given the difficulty in evaluating a Lorentzian vertex amplitude for a single choice of spins, it is a rather 
fortuitous result that all the spin dependence can be integrated out — indeed, as we shall see shortly, one 
obtains exact expressions in terms of the remaining and r e variables. 

We define the method explicitly as follows. Assuming that the original partition function is not affected 
by interchanging the order of integration, we can write: 

(2-6) Z A = f [J rdpf] I l[pjAf(pf)) (l[A e (pf) 

\f £ F Jo J \f£F J \e£E 



II II / 3 dX (e>») Y[ K Pf( X e+( 



(f,v), x e_(f,v) 



4 In this we follow the original Lorentzian model of Barrett and Crane and use exclusively the (0, p) representations from the 
principal series. While some proposals ]19l have been made which also include half integer labelled representations (fc, 0), most 
work has been with the original proposal. 

5 Some authors absorb the pj? factors from the measure into the definition of Af(pf); in the present work we keep the measure 
and face amplitude distinct. 
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n ( n / 3 das («.") ) ] (n f */) (in^/)) (u a m 

veV \e3v,e^e% J H *+ J J \f£F J ° J \f£F J \eeE / 

II ]J K Pf{ X e + (f,v),X e _ (/,«)) 

On the RHS of (2.6), all dependence on the pf variables is integrated out before performing the integrals 
over H+. As called for by the regularization, one integration with respect to H+ has been dropped from 
each vertex (the edge at v for which integration is dropped is denoted by eg) . To perform integration over 
the pf variables first, we consider the product of all the edge and vertex amplitude integrands, and organize 
this product into factors that each depend only upon a single pf variable. 

If we consider for the moment only the factors contributed by the vertex amplitudes A v given in (2.2), 
we note the following. Due to the form of the integrand of A v in terms of the functions (2.3), we see that 
for any given face / of the 2-complex all of the factors involving the spin pf can be grouped together into a 
vertex face factor Fj of the form 

v f sin(p f 4) sm(p f 4)..- sm(p f 4 Aogv{f) ) 
(2-7) F f (jp f ,<t>{) = , 

where the <\>[ denote distances on H\ and degv(/) denotes the number of vertices contained in the dual 
face /. With regard to the subscript i indexing the 4>{ variables, we assume that a numbering of vertices 
in every dual face / has been chosen; the i subscript runs over the values of the indexing map Nf(v) 
and the corresponding vertex v is found by inverting Nf(v), Explicitly, we choose for each / a bijective 
map of the form Nf : Vf — > {1, . . . , degv(/)}, where Vf is the set of vertices in /. Although we have 
suppressed the arguments of each <p{ in the formula above, for a given <f){ letting v — NJ 1 ^) we see that 

4>{ = $ ( x e- (f,v)i x e + (f,v) ) is a function of the variables x e _ff tV \ and x e+ tf :V y 

Although the face amplitude for the Perez- Rovelli model has no spin dependence (Af = 1), the edge 
amplitude (2.5) is a non-trivial function of the four spins variables associated to an edge of the 2-complex. 
While the eye diagram can be integrated to an exact expression in terms of the product of spin variables and 
the hyperbolic cotangent function [6J, the result does not have a form in which the spins pf appear in separate 
factors. However, we observe from (2.5) that before integration with respect to r e , the pf dependent part of 
the integrand is a product of factors of the form mn ^f re ' t Hence, if we collect factors from the different edge 
amplitude integrands that depend on the same pf, we have a product for each face which has the same form 
as the vertex face factor (2.7), but with the r e variables playing the role of the hyperbolic distance variables 
4>{ ■ Therefore, if the r e variables are integrated only after integrating with respect to the pf, we can define 
an overall face factor that takes into account contributions from both edge and vertex amplitude kernels. 
Our final form for the face factors are then 

(2.8) Ff(pf,<f>{, r e ) = ^L____ , 

p f 

where deg E (/) denotes the number of edges contained in the face / and the e(f,i) selects the ith edge 6 
contained in the face /. It should be noted that we have absorbed the pj terms from the measure into our 
definition of Ff by lowering the power of pf in the denominator by two. 

Having defined our face factors, we can now rewrite the partition function in terms of the integrals of the 
face factors Ff(pf,(p{,r e ) and products of sinh functions that depend only upon the r e and <fi{ : 




(2.9) Z A = ml dr. n II L, dx 




As with the face- vertex numbering Nf(v), we assume that for every face a numbering of its edges Af (/, e) has been chosen. 
We use an edge valued subscript for the r e as there is only one such variable for every edge. 
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Formula (2.9) is our explicit face factored formulation for the partition function of the Perez-Rovelli model; 
in Section 3 we shall prove it is equal to the original model for all non-degenerate triangulations and certain 
degenerate triangulations. 

From (2.9), we see that face factoring requires interchanging the order of integration among the pf, the r e , 
and the hyperboloid integrals — all of which are improper. In order to justify the interchange, it is sufficient 
that at least one of the orderings of integration is absolutely integrable. This follows from the Tonelli-Hobson 
test, a corollary of the well-known Fubini theorem in the case of integrals on an unbounded domain; see for 
example [Q. In Section 3 below, we justify this interchange for a general class of triangulated 4-manifolds 
by establishing absolute integrability for the face factored formulation of the partition function. 

We consider next an illustration of the face-factoring method for the 2-complex dual to a particular 
triangulation of the 4-sphere S 4 . 

Example face factoring of a Barrett-Crane model with trivial edge amplitude. For simplicity, we 
consider in this section a modified Perez-Rovelli model with (2.5) replaced by a trivial edge amplitude so 
that 

A f = 1, A e = 1. 
The standard Barrett-Crane vertex amplitude 



X e+ (f,v) ? X e _ (f,v) ) 



is used; recall from (2.3) that the kernel functions are 

sm(p f (/>(x,y)) 



K p (x,y) 



Pl ~' yy pfsvah{(j){x,y))' 

Consider the triangulation of S A resulting from taking the boundary of the 5-simplex. In this triangulation, 
there are 6 four-simplices, 15 tetrahedra, and 20 triangles. We shall work with the dual 2-skeleton of this 
triangulation, in which there are 6 vertices, 15 edges, and 20 polygonal faces; the relevant face factors 
Ff(pf ,(/){) will be deduced from counting and symmetry arguments. Because each vertex is contained in ten 
faces, there will be ten K Pf (x, y) factors per vertex for a total of sixty K Pf {x, y) factors. In our present case, 
all faces have an associated face factor of the same form by symmetry. Since by definition there is a face 
factor for each face / and there are twenty faces for the entire 2-complex, we see that each face factor has 
the form 

(2.10) FtbtA) = sin (P/^) sin ^^) sin ^^) ; 

Pf 

where the three (f>{ are associated to the three vertices in /. Having found the face factor, we consider the 
integral that appears in the face-factored form of the partition function definition given in (2.9): 

(2.1D r F fiPf , 4) <% = r sia(Pf4) sHpf4) sin(j)f4) dPf 

Jo Jo Pf 

= ^(sign(-<^{ +<!>{ + <t>l) - sign(-</>{ +<t>{- (j>{) - sign(</>{ +<!>{ + <t>l) + sign(^{ + <$>{ - <$>{)). 

As suggested above, this integral is not difficult to evaluate to an exact closed form in terms of the distance 
variables <f){, We note here that given the difficulty of evaluating a single Lorentzian vertex amplitude, it 
may be surprising that integrating over all the spin variables can be done in such a way that the result is a 
product of simple functions of the variables. 

We again emphasize that to equate the face factored form of (2.9) with the original formulation of the 
partition function, interchanging the order of the improper integrals needs to be justified. We won't provide 
such a justification for this example case — however, in Section 3 we prove that the face factoring method is 
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valid for the actual case of interest, the Lorentzian Barrett-Crane model with the Perez- Rovelli face and edge 
amplitudes. Before doing so, we shall give an exact closed form expression for the face factor integrations 
that arise in the Perez-Rovelli model. 

Closed form for the Perez-Rovelli face factor integrations. In the Perez-Rovelli model, the most 
general face factor is given by (2.8), and the integrals required to eliminate dependence on pf for a given / 
are of the form 

. v n^)(n«)^, 

Clearly, this integral can be expressed using a product of sine functions as 




K vef J Vee/ 



{ I Yl sinc(p f <t> f Nf{v) ) J j JJ sinc(p/r e ) J pf dp f . 



Were it not for the factor of pj due to the measure, one would be integrating a product of sine functions 

with arguments scaled by the 4>i and r e variables. If this were the case, a rather remarkable closed form 
given by Borwein et. al. in [TJl could be directly applied. This formula, which will be useful for our actual 
case as well, we give here as follows. 

Without loss of generality, let an arbitrary face / in the 2-complex be chosen. We define a new set 
of variables di so that di = 4>{ +1 for < i < degy(f) and = f e (/,i+i-dag v (/)) f° r degv(/) < i < 
degv(/) + degE(/)- As the di are all non-negative, we can apply from JHj the result: 



Jo \ l=0 Pf J z * n - 7e{ _ ljl} „ 

where n = degv(/) + degE(/) — 1 and 7 = (71, • • • , 7«) E {— 1, 1}™. The b 1 and e 7 are defined by 



b 7 = d Q + ^2 Ikdk, e 7 = \\ Ik- 

k=l k=l 

Remarkably, it turns out that this formula enables us to handle the Perez-Rovelli case as well — we simply 
differentiate both sides of (2.12) twice with respect to (any) one of the di parameters. On the LHS the 
differentiation is passed under the integral sign; choosing do as our parameter to differentiate we have: 



Note that differentiating sin(p/d ) twice with respect to do gives the desired factor of pj multiplying the 
negative of sin(p/do). Thus (2.13) provides a closed form expression for all the face factor integrals that may 
appear in the Perez-Rovelli model. In practice, one can also evaluate these types of integrals using symbolic 
integration with software such as Maple or Mathematica. 

3. Finiteness Results 

The finiteness of the Lorentzian Barrett-Crane partition function with the Perez-Rovelli choice for face 
and edge amplitudes was an important finding established 7 in 14] for arbitrary non- degenerate 2-complexes. 
In this section, we prove that the partition function is absolutely integrable for a more general class of 
2-complexes. This proof both justifies the interchange of improper integrals required by face factoring and 

7 See also |15| for further detail. 
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at the same time proves finiteness of the model for all non-degenerate triangulations and a certain class 
of degenerate triangulations. We characterize this class of degenerate triangulations and show how it is 
sensitive to certain changes in the face amplitude. Note that the proof we give in this section applies to 
closed two-complexes. The result also holds for 2-complexes with boundary; we describe in Appendix B how 
the proof of the closed case is changed to account for the presence of boundaries. 

We start by introducing an absolute bound on the Perez-Rovelli partition function in its face-factored 
form (2.9) above. 



(3.1) 



\Za\ < 




n n 



<j H 3 



dx 



(e,v) 




7rsinh (r e ) 



n n 



n 



v \f3v Sinh (<l>N f (v)) ) J \ f J ° 

Introducing the inequality (A. 5) proven in Appendix A: 

/ /-oo \ ( a f dc Sv(f) dcg E (/) 

( 3 - 2 ) n(/ o F (PfA{,r e ) d Pf j <n u i n 4 n 

Hence, if the RHS of the bound 



F f(Pf,<f>i,r e ) dp f 



~ dogv (/)+d<=g E 



r e(f,j) 



J 



(3.3) 



\za\< n 




n n 



dx 



(e,v) 



n(n sinh( ^ /w) 




/ deg t ,(/) deg E (/) 



n? n n 



1=1 




is finite then our face factored form of the partition function (and hence the original form) is finite. 

To integrate out the r e variables in our bound, for each e we collect together all the r e dependence into a 
single term of the form 



(3.4) 



o sinh (r e ) 



where a e > is the overall power of r e appearing in (3.3) due to contributions from the four faces that 
contain e. 

We claim that for a e > 1, quantities of this form are always finite. To see this we divide the integration 
as follows: 



I 



o sinh (r e ) 



■ dr e . 



o sinh (r e ) 



■ dr e 



hyz sinh (r e ) 



■ dr K . 



Consider the first term. Although the integrand is unbounded at r e = 0, sinh(r e ) — > r e as r e — > 0, so near 
zero the integrand behaves as r^^ a "\ which is integrable for a e > 1, but will otherwise diverge. Because 
^g- < sinh 2 (a;) for x > H^-, we can bound the second term as 



a sinh (r e ) 



■ dr e < 16 



i: 



dr e 



which is finite, as the polynomial r"" is exponentially damped as r e — ► oo. Thus all factors involving r e can 
be integrated to give finite factors if for each edge e the inequality a e > 1 holds. We will return below to a 
discussion of this a e > 1 condition and the constraints it places on which degenerate triangulations can be 
shown finite. 
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Next we need to check that integration over the hyperboloid variables always yields finite factors. Upon 
integrating out the r e as described, we have reduced the bound (3.3) to: 
(3.5) 



|2 " sc (n( JL L — ))(n (n ^ J J n (T * 



cn n l^(-) n 




„ J H 3 , 1 \ f"" sinh(ei>{ r , s 




1- 



do* v (/)-HU« E </> 



' do gv ,(/) + dcg E (/) 



where the r e integrations and all other constants have been absorbed into an overall constant C. 
From this expression, we see that to complete our proof we require 



(3.6) J] / dx 



I ™J(e,v) I n . t 

, vJh 3 / T-T smh f,, , J 

for each v to be finite. Note that these integrals are of the same form of the vertex amplitude (2), but 
with modified kernel functions. In recent work by Dan Christensen Q2]> a proof is given in which a class 
of integrals that includes those of the form (3.6) are shown to be finite. From this we can conclude that all 
hyperboloid integrations produce finite factors. 

Therefore, absolute integrability for our bound on the face factored formulation depends only upon finite- 
ness of the r e variable integrations, for which we found above that a e > 1 is a necessary and sufficient 
condition. 

We analyze the a e > 1 condition as follows. For all non-degenerate triangulations, degy(/) + deg B (/) > 6, 
hence the power of r e appearing in each face factor is at least (1 — |) = |. Since each r e is associated to 
a tetrahedron, it will appear in 4 face factors. Hence a e > (4 • |) = 2 > 1, and so we conclude that all 
non- degenerate triangulations are absolutely integrable. This justifies the face factoring method for these 
triangulations and gives an alternative finiteness proof to that of |14j . 

More generally, the condition on a e in terms of the vertex and edge degrees degy(/) and deg E (f) can be 

given as a e — J2f3e fl ~ dcg ^^ deg > 1. If we consider degenerate triangulations where deg B (/) > 2 
and deg E (f) > 2, this condition is equivalent to the following statement: 

Theorem 3.1. The Perez-Rovelli partition function Za is absolutely integrable if for every e £ A, deg v (f) + 
deg E (f) > 4 for at least one / containing e. 

This is a generalization of the finiteness proof of pi], which was limited to non-degenerate triangulations. 
We shall see next that some degenerate triangulations satisfy this condition and some do not. 

Finiteness for Degenerate Triangulations. In the case where a dual 2-complex contains a face / where 
degy(/) or deg E (f) is less than 3, it is dual to a degenerate triangulation. Such a "triangulation" is not 
strictly speaking a triangulation in the usual sense of a simplicial complex, as it contains 4-simplices whose 
intersection contains more than one tetrahedra. While visualizing 4-dimensional geometry is challenging, we 
can draw lower-dimensional analogs of degenerate triangulations; two examples are shown in Figure 3.1. 

Using Theorem 3.1, the 4-dimensional analog 8 of Figure 3.1(a) can be shown finite. Observe that although 
the interior degenerate edges (indicated with thickened lines) have a face where deg E (/) = 2 and deg y (/) = 2, 
they also have 3 faces where deg E (f) > 2 and degy(/) > 2, and so applying our criterion above we see that 
Za is finite. 

The four dimensional analog of Figure 3.1(b), in which two 4-simplices intersect along all five of their 
boundary tetrahedra, results in deg E (/) = 2 and deg y (/) = 2 for all faces inside the tetrahedron of the 
intersection; for edges dual to these tetrahedra a e = 4 • \ — 1 so integration of our face factor bound in the 
r e variable diverges. Hence our proof fails to show that such 2-complexes are finite. We remark here that 
our Theorem 3.1 is a sufficient condition but we have not shown it to be necessary — our work does not 
rule out finiteness of the face factored form for this type of triangulation. However, even if the face factored 




icgy (/) + dcg E i f ) 



In the analog, two four-simplices are glued along two of their tetrahedra. 
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form can be shown not to be absolutely integrable, the original form of the partition function may still be 
finite but only conditionally convergent. Currently, numerical computations are underway by the author to 
explore some of these possibilities for 2-complexes arising from these degenerate triangulations. 

A final observation relates to how possible changes in the face amplitude affect the integrability of de- 
generate triangulations. While the edge and vertex amplitudes we have used here are fairly well accepted 
and can be given some physical motivation, it may be of interest to consider alternatives to the Af(pf) = 1 
choice. If we allow ourselves to consider a face amplitude of the form 

MPf) = 4 
Pf 

for 7 > 0, we can construct a new bound (following the steps of Appendix A) in which the convergence 

condition becomes a e — J2f3e y- ~ dcg v (/Hdcg E (/) ) > ^' ^ e ana ^°S of Theorem 3.1 is than the requirement 
that 

4 

degv(/) + deg E (/) + -7 > 4 

for some / containing e. 

It is interesting to observe that for 7 > 0, this criterion is still met even if degy(/) = deg E (/) = 2 for 
all the faces of a tetrahedron, the case which fails for the original amplitude A/(pf) — 1 corresponding to 
7 = 0. Hence, with the original model such 2-complexes give a bound that is divergent but at the "margin" of 
integrability — any non-zero 7 results in a model where all degenerate triangulations satisfying deg E (/) > 2 
and deg E (f) > 2 are finite. 



4. Applications 

Generalization to the Livine-Oriti Causal Barrett-Crane Model. In the Livine-Oriti causal model 
developed in the kernel functions for the vertex amplitude A v are modified to take into account face 
orientations s{, <E { — 1,1} relative to a vertex v that arise when imposing causality. The causal kernel 
replacing the original kernel (2.3) is given as 

(4.1) K Pf (x,y) = ; . 

p f smh(<t> T Nfiv) (x,y)) 

From the form of this kernel, we can write the vertex face factors as 

(4-2) ^(P/^) = T^7TlI^' W - 

Pf vef 

The integral / °° Fj{pfi4>{) dpf diverges in the general case, because the numerator approaches a constant 
as pf — > while the denominator approaches zero at least linearly; this gives a divergent integral for 
degv(/) > 1, which applies to any triangulation. Hence, contributions to the overall face factor from face 
and edge amplitudes must be relied upon to give a finite face factored partition function. 
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Convergence criterion for the face amplitude: In the edge amplitude remains the eye diagram function 
when passing to the causal case. Taking the causal vertex and edge amplitudes as fixed, we show in this 
section how the face amplitude has to be carefully tuned to avoid divergence of the face factored formulation. 
Let us consider the full causal face factor with the edge amplitude kernels and a face amplitude of the form 
Af(pf) — Pj for a real number (3: 

(4.3) f>(p/,tf,r.) = dcmU)+ i m[i) ^ n^'^ W ( Il sin (^) ) • 

Pf \vef j \eef ) 

Observing that the product of degE(f) sine functions can cancel against the p^ CSE ^- ) at pf = 0, to avoid 
divergence we arrive at the condition 

deg v (/) - 3 < p. 

In addition, if (3 > degv(/) + degE(f) — 2 then Ff(pf, (j>{ ,r e ) for general <f> fails to have a limit as pf — > oo, 
and so its integral is not defined. Hence a necessary finiteness condition for a given f3 is 

(4.4) deg v (/) - K + 2 < deg v (/) + deg E (f). 

For = as has been proposed, (4.4) becomes degv(/)< 2 and degv(/) + degE(f) > 2; thus we see that 
the former condition restricts us to degenerate triangulations. For j3 > 1 however, degenerate triangulations 
are completely ruled out. Note here that in (4.4) the parameter (3 is bounded from both above and below. 

We stress here that failure of convergence for the face factored partition function does not imply divergence 
of the original form; it may still be that the original partition function is conditionally integrable, precluding 
the identification of the two forms. 

To summarize, we find that for the Oriti-Livine model there are some rather strong constraints on Af(p/) 
necessary for the existence of a finite face factored formulation. An interesting goal for future work is to de- 
termine whether the model with original integration ordering escapes these constraints by being conditionally 
convergent, or if it in fact is divergent in the same cases as are divergent in the face factored formulation. 

Numerical Applications. Within the context of Monte Carlo methods applied directly to the partition 
function, the face factored formulation is advantageous in several respects. The integrations eliminated 
involved a highly oscillatory product of sine functions; the integration of these functions leads to much 
smoother functions of the remaining variables. In fact, it can be shown that the remaining functions are 
always non-negative [H]. The features of smoothing and positivity are both desirable in Monte Carlo 
methods for improving the accuracy achievable for a given number of samples. Currently, computations are 
underway to evaluate the partition function for sufficiently small 2-complexes where the number of integration 
dimensions is still not prohibitive. 

Apart from direct Monte Carlo methods which have severe scaling limitations, the face factored form of 
the partition function may be amenable to statistical mechanical methods similar to those used in lattice 
quantum field theory computations. In such computations, one avoids calculating the partition function 
itself, but rather generates a sequence of samples in configuration space which has the same distribution as 
that determined by the partition function. In this way, expectation values of observables can be found to 
good approximation without calculating the entire partition function. The Metropolis algorithm is one such 
method that was successfully applied t 3.; in calculating expectation values for the Riemannian Barrett-Crane 
model. 

Many statistical mechanical methods require that the partition function amplitude 9 be non-negative for all 
configurations (i.e. the partition function amplitude is the exponential of some real- valued action). Clearly, 
before face factoring the Perez-Rovelli amplitude as a function of all the spin variables fails to have such 
positivity due to its definition in terms of sine functions. However, after the pf variables have been integrated, 
it is rather remarkable that for closed spin foams the amplitude as a function of tetrahedral variables is always 
non-negative JH- The potential for applying lattice field theory techniques to spin foam models using "dual" 
tetrahedral variables has been emphasized by Pfeiffer in L 2*Dj, and may open the door to computing with 
much larger 2-complexes. 

9 The integrand of the partition function a function of pf variables in the original formulation and a function of if? and 
r e variables associated to tetrahedra in the face factored formulation. 
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On a cautionary note, the physical interpretation of a model defined purely in terms of tetrahedral variables 
seems ambiguous at this time. Based upon connections to the Regge action in certain limits and general 
geometric considerations, it has been argued 10 that the hyperboloid variables can be interpreted as the time- 
like normals to the tetrahedra (which are purely spacelike by the Barrett-Crane construction). While this 
viewpoint is certainly promising, it is unclear to the author how these interpretations can be related back to 
the canonical approach. 

5. Conclusions 

In this paper we have shown how the Lorentzian Barrett-Crane spin foam model can be reformulated 
so that integration with respect to the spin variables p/ is performed before all other integrations to yield 
exact expressions in terms of tetrahedral variables. Absolute integrability was proven for the face factored 
formulation, justifying the interchange of integration and establishing finiteness for a large class of triangula- 
tions. This class includes all non-degenerate triangulations and a certain class of degenerate triangulations, 
extending the finiteness proof of [2j. In addition to the face-factoring transformation, an essential ingredient 
of this proof is a recent finiteness result due to Christensen ,12| for a large class of lOj-like integrals. 

We have also described how the type of degenerate triangulations proven finite by our method depends 
upon the choice of face amplitude. Given this result, it may be worthwhile to investigate more thoroughly 
the physical significance of the degenerate triangulations. While the method can also be applied to the causal 
Livine-Oriti model, we found that given their choice of edge and vertex amplitude, there is a rather stringent 
constraint between the form of the face amplitude and the types of triangulations that allow a finite face 
factoring formulation. 

An interesting feature of the Barrett-Crane model revealed by our work is that edge and vertex amplitudes 
contribute kernels to the face factors of the same form; this is no accident and can be traced back to the 
diagrammatics used to derive the model [6J. Hence, an obvious question is whether or not generalizations 
or modifications of the Barrett-Crane model may still allow a face factoring formulation. Specifically, one 
could consider whether the Barrett-Crane type model proposed in [Tflj . with mixed representations, may be 
susceptible to a face factoring treatment. Moreover, it may be interesting to investigate whether any analog 
of face factoring is possible for spin foam models not of the Barrett-Crane type. 

Acknowledgement. The author would like to thank Dan Christensen, Igor Khavkine and Joshua Willis for 
very useful discussions and suggestions. The author also acknowledges the financial support of NSERC and 
OGS. 

Appendix A. Product Bound for General Face Factor Integrations 

In this appendix we find functions of r e and (f>{ that bound from above the possible integrated face 
factors that arise in the Perez-Rovelli model. The bounding functions have the simple form of a product 
of r e and <f>{ variables raised to a power that depends only upon the number of edges and vertices that 
the face contains. These bounding functions are used in Section 3 to show the absolute integrability of the 
face-factored formulation. 

To begin, we recall that the most general face factor that can arise for the Perez-Rovelli model has the 
form 



(A.l) 



3in(p/<ft{ ) ■ ■ ■ sin(p/<^ 08y(/) ) sinfr/T-e^i)) ■ ■ ■ sin(p / r e(/4cgE(/)) ) 

do gv .(/)+dog E (/)-2 



where the r e and (f>j are non-negative real numbers coming from edge and vertex amplitudes, respectively. 
To see this, note simply that in both the edge and vertex amplitude each sine function in the numerator is 
accompanied by a factor of pf in the denominator. The measure then contributes a factor oi pf, reducing 
the degree of the denominator by 2. Based on this form, we can immediately establish two types of bounds: 

^deg v (/) deg E (/) 



(A.2) Ff( Pf , 



< 



p f 4>{ ■ ■ ■Pf<pj CSvif) Pfr e (f,i) ■ ■ ■P/^(/,dcg E (/)) 

deg v (/)+deg B (/)-2 
P f 



n 



n r <f>j) 



pf 



10 



For instance in a dual variable framework is also developed in some detail for the Euclidean case in [2^ 
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and 



(A.3) 



F f(Pf,<l>i, r e) 



sm( Pf <f>{) ■ ■ ■ sm(p f (P f degv{f) )sm(p f r eUA) ) ■ ■ • sin( P/ r e(/;degE(/)) ) 



dcg v (/)+dcg E (/)-2 
P f 



< 



deg v (/)+dcg E (/)-2- 



The first inequality (A. 2) is valid for pf > 0, while the second bound (A.3) can only be used for pf > 0. 
Integrating both sides of these inequalities we can form the following estimate for the integrated face factor: 



,-oc i-M / dc Sv(/) dcg E (/) \ 

(A.4) / F f {p f ,4,r e ) d Pf < [] 4 II r «M )Pf d Pf + 

Jo Jo y i=1 j=1 j Ja 



u p deg v (/)+deg E (/)-2 



, /deg v (/) deg E (/) \ 

^ n 4 n 



AT 



1 



M dc Sv(/)+dcg E (/)-3' 



for any M > 0. Since we are integrating pf for a given set of the 4>i and r e , we can make this M a function 
M((f>{ , r e ) of these variables. We make the following choice: 



M{4,r e )= ( n ^ n 

i=l j=l 



which upon substituting into (A.4) gives: 



(A.5) 



„oo -, /deg v (/) dog E (/) 

/ Ff(pf, (f>{ , r e ) d Pf <-[ n 0f n 



1 



M d °Sv(/)+dog E (/)-3 



( /deg v (/) dcg E (/) \ 1 do gv ,(/)+deg E (/) /deg v (/) dcg E (/) \ 1 dog v (/) + deg E (/) 

U n ^ n + n ^ n ^ 

»=i j=i / V i=i i =i 



4 
3 



M<'!Sv(/) dcg E (/) \ ' r ■ ■' - 

II II r e(/,j) 

3 = 1 



Although we assume M > for this derivation to be valid, if any of the <j>j or r e are zero than the original 
integrand (A.l) and hence the integral over pf vanishes; hence the bound (A.5) is valid for any choice of <p{ 
and r e . As it relates to the discussion at the end of Section 3, we mention here that for a face amplitude of 
the form Af(pf) = the above proof can be generalized to arrive at a bound: 



/>oo / -, \ /deg v (/) dcg E (/) 



' do gv ,(/)+dog E (/) 



3) 



by following essentially the same steps. 
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Appendix B. Finiteness for spin foams with boundary 

As a path integral theory of quantum gravity, a spin foam model defines a transition amplitude from one 
3-geometry to another for a given "history" — a triangulated 4-manifold whose boundary is the union of the 
ingoing and outgoing 3- geometries. 

In the case where the 2-complex is dual to a triangulated 4-manifold with boundary, a straightforward 
extension of our finiteness proof for the closed case can be given. The spin variables associated to each face in 
the boundary are taken as boundary data, so these are not integrated in defining the partition function. As 
well, for dual edges and faces coming from boundary tetrahedra and triangles, one needs to take the square 
root of the edge and face amplitudes so that the partition function multiplies when two histories are "glued" 
along a common boundary. We shall denote the sets of edges and faces dual to the boundary tetrahedra 
and triangles as dE and dF, respectively. Rather than working directly with the square root of the edge 
amplitude ©4(pi, . . . ,f> 4 ), in this proof we use the following bound: 



(B.l) x/©4(pi, • • • ,Pi) < 1 + e 4 (pi, • • • ,Pi). 

Applying this bound to the partition function of a 2-complex A with boundary we have 
(B.2) 

Z A = j ■•■ f ( I] a/©4(pi,...,P4)) ( 11 ©4(pi,...,p 4 ) J ( HMpu---,Pw)) \{Pfdp } 

yO JO _ \„ c aij / \ / / Map 



fidF 




s \eedE / \efdE / \veV / f£dF 



4 (pi, • • • ,Pi) j Y[ 9 4 (pi, • • • ,Pa) j J| A„(pi, . . . ,pw) J J| Pf dpf. 

J \e(£dE ) \veV ) f£dF 



f$dF 



Considering the RHS of (B.2), we see that expanding the product of 1 + 84 functions results in a sum 
of partition functions in which boundary tetrahedra are either trivial (A e = 1) or have the eye diagram 
amplitude A e = 4 (pi, ■ ■ ■ ,Pa)- We will show that any partition function of this form is finite. 

First we handle the integrals with respect to the edge variables r e . For the given 2-complex A, let 
a partition function be chosen where the amplitude for any boundary edge is either A e — 1 or A e = 
4 (pi, . . . ,Pa), and all other amplitudes agree with the Perez-Rovelli model. We can extend the face factoring 
method to a partition functions of this form as follows. We begin by defining a boundary face factor Ff 
associated to any face on the boundary as: 

sin(p / ^{)sin(p/^)---sin(p/^ (/) ) sin(p/n) sin(p/r 2 ) • ■ • sin(p/r deg e (/) ) 

(B.3) F f (p f , 0. , r e ) = de gv .(/) + dc g e (/) 

p f 

where the function deg®(/) counts the boundary edges of type 04 contained in /. Observe that since the 
Pf are fixed on the boundary rather than integrated over, the measure factor of pj does not enter into the 

definition of Ff(pf, </>{ , r e ); the functions Ff(pf, 4>j ,r e ) and not their integrals with respect to pf multiply 
into the partition function directly. The analog of the inequality (A. 2) for faces in the boundary can be given 

as 

jp8, jJ ^ * PftiPfd ■ ■ ■Pf ( i ) L v U)Pf ri Pf r2 ■ ■ -P/ r degg(/) ( de T\ f \f de f\ f) 
(B.4) F f (pf,<p{,r e ) < = Yl Yl re{fJ) 

Pf \ »=i i=i 

In the open spin foam case, the inequality (3.2) is modified to include the boundary face factors: 

/*°° 1 
IJ / Ff{p f ,4,r e ) d Pf \F^{pf,4>lr e ) 

f£dA J ° fedA 

i-dcg^+dcg^/) /deg v (f) degg(/) 

11 11 4 n m«) 

fedA \ i=l j=l 
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For tetrahedra in the interior, finiteness goes through as in the closed case. For tetrahedra on the boundary, 
there is an integration over r e for any edge e of type 84; for such edges one collects all the r e dependence to 
give factors of the form: 



rt 

■ dr e 



sinh 2 (Y e ) 

As this integral is clearly finite, integrating with respect to the r e for all e in the interior and boundary 
always yields finite factors. For boundary tetrahedra where Ae = 1, there is no r e variable to be integrated 
so finiteness of the bound is not affected. 

The final step is to show that the result of integrating the partition function bound over the hyperboloid 
variables is always finite. For a 4-simplex which has all of its faces on the interior, we recall from the proof 
of the closed case that integration with respect to hyperboloid variables factors into products of the form 




' dcg v (f) + .h-± E if) 



f3v sinh(^ /W ) 

each of which can shown finite by invoking the result of 22 • F° r 4-simplices that have one or more faces in 
the boundary, one instead has 

n / 1 n [ ,iv)) f n t . 

which again are shown to be finite in ^2] — thus we have checked all hyperboloid integrations lead to finite 
factors. 

As our bound (B.2) is a finite sum of partition functions which are finite, we conclude that the Perez- Rovelli 
Lorentzian partition function is finite for all 2-complexes dual to triangulated 4-manifolds with boundary, 
on the condition that the interior faces meet the same criteria with respect to degy(/) and deg B (/) as were 
found for the closed case. 
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